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Abstract 

Background: CMR allows investigating cardiac contraction, rotation and torsion non-invasively by the use of 
tagging sequences. Three-dimensional tagging has been proposed to cover the whole-heart but data acquisition 
requires three consecutive breath holds and hence demands considerable patient cooperation. In this study we 
have implemented and studied k-t undersampled cine 3D tagging in conjunction with k-t PCA reconstruction to 
potentially permit for single breath-hold acquisitions. 

Methods: The performance of undersampled cine 3D tagging was investigated using computer simulations and 
in-vivo measurements in 8 healthy subjects and 5 patients with myocardial infarction. Fully sampled data was 
obtained and compared to retrospectively and prospectively undersampled acquisitions. Fully sampled data was 
acquired in three consecutive breath holds. Prospectively undersampled data was obtained within a single breath 
hold. Based on harmonic phase (HARP) analysis, circumferential shortening, rotation and torsion were compared 
between fully sampled and undersampled data using Bland-Altman and linear regression analysis. 

Results: In computer simulations, the error for circumferential shortening was 2.8 ±2.3% and 2.7 ±2.1% for 
undersampling rates of R = 3 and 4 respectively. Errors in ventricular rotation were 2.5 ± 1.9% and 3.0 ±2.2% for 
R = 3 and 4. Comparison of results from fully sampled in-vivo data acquired with prospectively undersampled 
acquisitions showed a mean difference in circumferential shortening of -0.14 ± 5.18% and 0.71 ±6.16% for R = 3 
and 4. The mean differences in rotation were 0.44 ± 1.8° and 0.73 ± 1.67° for R = 3 and 4, respectively. In patients 
peak, circumferential shortening was significantly reduced (p< 0.002 for all patients) in regions with late 
gadolinium enhancement. 

Conclusion: Undersampled cine 3D tagging enables significant reduction in scan time of whole-heart tagging and 
facilitates quantification of shortening, rotation and torsion of the left ventricle without adding significant errors 
compared to previous 3D tagging approaches. 

Keywords: Motion analysis. Myocardial tagging, k-t undersampling. Principal component analysis. 
Image reconstruction 



Background dyssynchrony [2]. Among the various mechanical 

Patients with myocardial infarction (MI), coronary artery parameters, cardiac strain and torsion have been shown 

disease (CAD) or impaired signal conduction in the valuable clinical parameters in patients with MI [3] or 

heart such as left bundle branch block (LBBB) suft'er aortic stenosis [4,5]. Cardiovascular magnetic resonance 

from reduced cardiac function. In patients undergoing (CMR) offers a functional imaging modality to assess 

cardiac resynchronization therapy (CRT) [1] it is of great cardiac motion pattern and synchrony of contraction 

interest to generate a mechanical map indicating local [6,7] non-invasively. 

Up to date multiple approaches have been used to 
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encoding [10-12], tagging by spatial modulation of 
magnetization (SPAMM) [13] and complementary 
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spatial modulation of magnetization (CSPAMM) [14]. 
These methods can either be applied in conjunction with 
two dimensional (2D) or three dimensional (3D) im- 
aging. In 2D acquisitions, multiple slices are imaged 
along the left ventricle (LV) and strain maps are calcu- 
lated slice-by-slice. Two-dimensional acquisitions re- 
quire additional techniques in order to compensate for 
through-plane motion. To this end, slice following 
[15,16], acquisition of additional orthogonal slices [17] 
or the encoding of through-plane displacement [18,19] 
(zHARP) have been proposed. While reconstruction of 
3D strain patterns from two bi-planar acquisitions does 
require interpolation, slice-following techniques provide 
only a projection of true 3D motion onto a 2D subspace. 
Using zHARP additional gradients in through-slice dir- 
ection are applied to estimate through-plane displace- 
ment by solving a set of linear equations [18]. 

To circumvent the need for interpolation and slice- 
following techniques, true 3D tagging has been proposed 
[20]. By applying modulation of magnetization in all 
three spatial directions, through-plane motion is intrin- 
sically captured and hence no slice-following is needed. 
Furthermore, 3D acquisition provides full LV coverage 
and yields an intrinsically higher signal-to-noise (SNR) 
compared to 2D imaging. While cine 2D acquisitions 
easily fit into a single breath hold, 3D acquisitions re- 
quire multiple breath holds with long duration [20]. Fur- 
ther approaches attempted to shorten scan time to four 
breath holds of 21 R-R intervals each [21] by only sam- 
pling around the harmonic peaks in k-space. These 
efforts, however, led to a lower spatial resolution in the 
two phase encoding directions. A different approach was 
proposed by Zhong et al. [22] using a 3D DENSE acqui- 
sition covering the entire left ventricle during free 
breathing. However, acquisition durations of up to 
20 min [23] depending on navigator efficiency were 
reported and hence feasibility for routine clinical use 
may be questioned. Rutz et al. [24] introduced an accel- 
erated 3D CSPAMM method only requiring three navi- 
gator gated breath holds of 18 R-R intervals each. This 
implementation was tested in a clinical study with 
patients suffering from MI and LBBB [25] providing 3D 
maps of synchrony and magnitude of contraction. To 
address issues of different breath hold levels the data 
obtained from respiratory navigators was used to correct 
for potential offsets. In their study the respiratory navi- 
gator was placed onto the right hemi-diaphragm and 
thereby respiratory induced displacement of the heart 
was approximated according to a linear relationship be- 
tween liver and heart displacement [26]. The accuracy of 
the linear translation from motion of the lung liver inter- 
face to the position of the heart is still being debated. 
Nehrke et al. [27] reported a strong correlation of the 
displacement of the right hemi-diaphragm and the heart. 



but also found significant subject dependent variability 
in the correction coefficients especially comparing 
inspiration and expiration. Subject dependency of 
the translation of breathing induced liver motion to 
bulk motion of the heart has been confirmed by 
Moghari et al. [28]. 

In general, data acquisition can be accelerated by 
undersampling in spatial and temporal dimensions. 
Among the various approaches, two strategies have 
gained significant attention. Compressed sensing (CS) 
[29] employs non-linear reconstruction methods to re- 
cover information from randomly or pseudo-randomly 
undersampled data. Inherently, compressed sensing 
algorithms require incoherent sampling and hence be- 
come applicable if a sufficient number of phase-encodes 
exists. In case of one-dimensional tagging preparation, 
the number of Cartesian phase-encodes orthogonal to 
the tag direction can be greatly reduced providing an ef- 
ficient and simple way for scan time reduction. In conse- 
quence the degrees of freedom to generate random 
sampling patterns become very limited and the applica- 
tion of compressed sensing appears less favourable in 
this particular application. 

The second acceleration strategy involves uniform spa- 
tiotemporal undersampling in conjunction with linear 
reconstruction algorithms. In k-t BLAST and k-t SENSE 
[30] low spatial but full temporal resolution training data 
is used to unfold signal aliasing resulting from data 
undersampling. The drawback of these methods relates 
to temporal filtering if undersampling rates increase. To 
address this issue. Principal Component Analysis (PCA) 
of the spatial-temporal frequency domain data was 
introduced and results obtained with k-t PCA show 
improved temporal fidelity [31]. 

Although k-t undersampling has been extensively ap- 
plied in CMR including cine and real-time imaging 
[32,33], perfusion [34-39] and phase contrast imaging 
[40-42], only few attempts of applying k-t undersampling 
to tagging have been reported [43,44]. 

The objective of the present study was to implement 
and test k-t undersampled whole-heart 3D CSPAMM 
tagging for rapid assessment of cardiac motion. The per- 
formance is demonstrated on simulated data, data 
obtained in healthy subjects and in patients with myo- 
cardial infarction. 

Methods 

k-t PCA 

To reconstruct undersampled 3D tagging data, k-t PCA 
[31] is used. Similar to k-t BLAST [30], acquired data is 
divided into 1) training data with low spatial resolution 
(in phase encoding directions ky / k^) but full temporal 
resolution and 2) k-t undersampled data with high 
spatial resolution and full temporal resolution. The 
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training data Ptrain(^>t) is Fourier transformed to be 
represented in the spatial-temporal frequency domain 
{x-f). Using principal component analysis PtraiJ^J) is 
then decomposed into a basis of temporally dependent 
functions b(fj) corresponding to the principal compo- 
nents (pc) and spatially dependent weighting coefficients 
^trainM 1^ x-pc space accordiug to: 



Ptn 



^{Xijj) = b(fj)Wtrain{Xif 



(1) 



The aliased signal at point (x, fj) resulting from R-fold 
undersampling can be written as: 



P alias {xjj) =^b{fj) "^{Xif 



(2) 



with w(Xi) denoting the spatial weighting coefficients of 
the unaliased image. Hence the aliased signals Paiias(^>f) 
can be expressed as: 



- alias 



P alias i^^fl) 
Palias{xj2) 

Paliasi^jN) 



Ewx with Wx 



w'^{xi) 

w'^{X2) 


and E = 




w'^{Xr)_ 







(3) 



Finally, the unaliased spatial weighting coefficients are 
obtained by solving: 



(4) 



where 0 represents an estimate of from training data, 
W denotes noise variance, H conjugate transpose and + 
the Moore- Penrose pseudo-inverse. 

Computer simulation 

All simulations were performed in Matlab (The Math- 
Works, Natick, MA, USA). Three orthogonal stacks with 
line tagging modulation in readout direction were gener- 
ated (Figure la). The CSPAMM method [14] was simu- 
lated to avoid tag line fading. The model consisted of a 
contracting left ventricle as well as static tissue repre- 
senting chest wall and liver. Circumferential shortening 
and rotation as measured in a healthy subject at basal 
and apical level was linearly interpolated along the long- 
axis to create three-dimensional motion data. Peak cir- 
cumferential shortening was 18.8% and 17.8% for base 
and apex, respectively. Peak rotation was -3.2° and 10.3° 
for base and apex (Figure Ib-d). Longitudinal shortening 
obtained from the same in-vivo subject was incorpo- 
rated. Simulations of undersampled data acquisition 
were compared to fully sampled simulated data sets with 
equivalent spatial and temporal resolution. The matrix 
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Figure 1 Numerical simulation. Three stacks with orthogonal line 
tag pattern were simulated (a). As input for the model, longitudinal 
shortening (b), rotation (c) and circumferential shortening (d) 
obtained from one in-vivo acquisition were used. 



size was set according to practical values [24,25] 
(Table 1). Gaussian noise was added to k-space data be- 
fore undersampling resulting in a SNR of 25 prior to 
undersampling. Both undersampled and training data 
were extracted from the computer model (Figure 2b). 
Undersampling rates of R = 3, 4, 5 and 8 were simulated. 
In all simulations, five training profiles were used in ky 
and kz direction resulting in a total of 25 training 
profiles. 

In order to study regional wall motion abnormalities, 
myocardial infarctions of different severity were simu- 
lated. To this end, radial shortening in the lateral sector 
was changed from 100% (no infarction) to 0% (com- 
pletely static) in steps of 10% (Figure 3). The reduction 
of myocardial motion in the infarcted zone was applied 
transmurally along the entire long axis of the left ven- 
tricle. Reduction of radial shortening directly reduced 
the circumferential contraction. In order to assure a 
smooth transition between infarcted and healthy tissue, 
the infarcted tissue was continuously "attached" to the 
adjacent healthy tissue, by reducing the motion damping 
factor continuously over a sector of 40° on both ends of 
the infarcted region. 

In-vivo measurements 

Eight healthy subjects (7 male, age: 27.5 ± 3.5 years) and 
five patients with myocardial infarction (5 male, age: 
54.8 ± 5.9 years, 1 patient with acute myocardial infarc- 
tion post percutaneous coronary intervention (PCI), 4 
patients with chronic myocardial infarction) were stud- 
ied. Imaging was performed on a 1.5 T Philips Achieva 
System (Philips Healthcare, Best, the Netherlands) using 
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Table 1 Parameters used for numerical simulation 



Undersampling factor 


Matrix size kx x ky 


Slices kz 


Heart phases 


Training profiles (kyXkJ 


Rnet 
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24 
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15 


24 
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65x16 


15 


24 


25 (5x5) 


44 



a 5 channel cardiac receiver array. Written informed 
consent was obtained from all subjects and the protocol 
was approved by the institutional review and ethics 
boards. 

Tagged images were acquired using an ECG triggered 
multi-shot EPI sequence with sequence parameters listed 
in Table 2. As illustrated in the sequence diagram in 
Figure 2 the tagging preparation was applied within 
10 ms after the detection of the R-wave. As the first RF 
pulse used for tagging preparation was applied select- 
ively in phase encoding direction, signal outside the 
field-of-view was suppressed after CSPAMM subtraction 
and hence reduced field-of-view acquisition could be 
performed [16]. The imaging sequence was repeated to 
acquire three orthogonally tagged stacks. For each stack 
21-24 heart phases were recorded. 

Fully sampled data was acquired in three consecutive 
breath holds. To guarantee similar breath hold position 
for each tagged stack a pencil beam respiratory navigator 
placed on the right hemi-diaphragm was used. The pos- 
ition of the lung liver interface was monitored at the be- 
ginning of each breath hold and displacements within a 



5 mm gating window were accepted. The breath hold 
was repeated if the subjects breath hold level was not 
within the given window. Retrospective undersampling 
was applied with reduction factors of R = 3 and 4. Pro- 
spectively undersampled data were acquired according to 
the schematic shown in Figure 2b in an additional breath 
hold. To facilitate comparison, all imaging parameters 
except for undersampling factors were kept identical for 
fully sampled and undersampled data acquisitions. For 
undersampling factors R = 3 and 4 the corresponding sets 
of parameters are given in Table 2. Prior to image recon- 
struction, the 5-channel coil array data were compressed 
into a single virtual coil data set using the array compres- 
sion method proposed by Buehrer et al. [45]. 

Data analysis 

Epicardial and endocardial contours were manually 
drawn for each slice of interest and a midmural contour 
was calculated. This contour was tracked over time 
using the peak combination harmonic phase (HARP) 
[46,47] algorithm implemented in an in-house software. 
Rotational motion and circumferential shortening was 
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Figure 2 Sequence diagram (a) and sampling pattern for k-t 
undersampling (b). After detection of the R-wave the tagging 
preparation is applied, followed by two repetitions of undersampled 
data and training data acquisitions for each stack. 
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Figure 3 Comparison of peak circumferential shortening from 
undersampled (R = 3) and fully sampled simulated data. 

Myocardium was divided into three groups: infarction (I), which 
covers a sector of 80° (22% of myocardium), healthy tissue (H) and 
infarct adjacent (A) tissue "attaching" infarcted to healthy tissue over 
a sector of 40° (1 1% of myocardium). For infarcted tissue, different 
degrees of immobilized tissue ranging from 0% to 100% were 
simulated. The dashed line indicates the reference diagonal with 
slope 1. 
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Table 2 Imaging parameters for the acquisition of fully 



sampled reference and undersampled data (R = 3 and 4) 



1 InHorcamnlinn fa^tof 
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IVICILI lA OlZ-C l\x l\y 


28 X 1 5 
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Phase-encodes 


15 


15 


Receiver bandwidth 


314 Hz/pixel 


314 Hz/pixel 


Spatial resolution 


3.8x7.2x7.2 mm^ 


3.8x5.4x7.2 mm^ 


Heart phases 


21 


24 


Temporal resolution 


30.3 ms 


30.3 ms 


EPI factor 
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TE^R 
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Rnet 
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Tag line distance 
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calculated from tracked contours as previously described 
by Ryf et al [48]. Curves of circumferential shortening 
and rotation over the cardiac cycle were fitted by 8^^ 
order polynomials, and maxima were found by estimat- 
ing roots of the derivative. From the mean rotation over 
the contour, torsion was derived by taking the difference 
in peak rotation at the most apical and the most basal 
level [49]. Analysis of peak circumferential shortening 
and rotation was performed in six sectors per slice in 9 
slices for computer simulation and 8-9 slices in in-vivo 
experiments. Comparisons were performed sector-wise 
and are reported for the entire LV. 

For simulation the tracking results from fully sampled 
data were used as ground truth reference. For in-vivo 
measurements the mean of paired data points was used 
as reference. In order to facilitate comparison the initial 
contour was kept the same for the fully sampled refer- 
ence and the retrospectively undersampled data in simu- 
lation and in-vivo. In the comparison of fully sampled 
data with prospectively undersampled acquisitions, con- 
tours were redrawn to account for changes in breath 
hold position and/or patient motion in-between scans. 

Analysis of simulated data was done by calculating 
the relative difference defined as: error rei = 
\max{motionii=i ) — max{motioniiyi) \ / \ max{motionR=i ) — 
min{motionR=i)\ , The fully sampled data was used as 
reference. Relative differences ± one standard deviation 
are reported in % for circumferential shortening and 
rotation. In order to estimate the correlation between 
undersampled and reference data linear regression was 
performed correlating circumferential shortening and 
rotation of reference data and undersampled data. Re- 
gression slope, offset and the corresponding 95% 



confidence interval were estimated as well as the correl- 
ation coefficient R^ and the standard error of the esti- 
mate (SEE). The SEE is given in % of the range of 

motion and is defined as SEE = ^J ^^^-i ~ 
with N being the number of points, Xi denoting meas- 
urement data and x{ corresponding values derived from 
linear regression. For in-vivo imaging, peak circumferen- 
tial shortening and rotation were compared using Bland 
Altman analysis. Mean differences and the 95% levels of 
agreement corresponding to 2 standard deviations (2SD) 
are reported. Linear regression was performed on peak 
circumferential shortening and rotation. Regression 
slopes, offsets, correlation coefficients and the SEE are 
reported. 

Comparing reference data and retrospectively under- 
sampled data, the HARP tracking performance of con- 
tour coordinates was studied based on contour vertex 
definitions in polar coordinates. Similar to the analysis 
of peak circumferential shortening and rotation. Bland 
Altman and linear regression analysis were performed. 

Patient data are reported using bulls-eye plots of peak 
circumferential shortening, and late gadolinium en- 
hancement (LGE). The centre of the bulls-eye plots 
represents the apex and the outer ring the base of the 
LV. For LGE images the average myocardial signal inten- 
sity per slice and sector normalized to the signal inten- 
sity measured in infarcted tissue is shown. Sectors were 
grouped into two groups: 1) non-viable sectors having 
more than 50% of the area presenting LGE and 2) viable 
sectors having less or equal to 50% of the area present- 
ing LGE. Peak circumferential shortening was estimated 
for all sectors in each group and the median, 50^^ per- 
centiles and 90*^ percentiles were estimated and pre- 
sented in box-plots. Statistical significance of the 
differences in peak circumferential shortening between 
both groups were estimated by a two-tailed Wilcoxon 
signed-rank test. The results were Bonferroni corrected 
for repeated testing. A p-value less than 0.05 was consid- 
ered statistically significant. 

Results 

Computer simulation 

Normalized differences in peak circumferential shorten- 
ing and rotation between undersampling factors of R = 3, 
4, 5 and 8 and fully sampled data used as reference are 
shown in Figure 4a-b. Errors in circumferential shorten- 
ing and rotation were below 5% for undersampling fac- 
tors of R = 3-5 and R = 3, 4, respectively. At R = 8 errors 
in circumferential shortening and rotation were 2.5 and 
3.1 times greater compared to data obtained with R = 4. 
The differences in torsion between undersampled and 
fully sampled data were -0.13°, 0.42°, 0.95° and 1.85° for 
R = 3, 4, 5 and 8, respectively. 
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Figure 4c-j shows regression slope, offset, correlation 
coefficient and SEE for R = 3, 4, 5 and 8. The 95% confi- 
dence intervals are presented as error bars for fitted 
slopes and offsets. 

Figure 3 shows the result of motion tracking in the 
presence of simulated infarction. Peak circumferential 
shortening in an equatorial slice is plotted for the refer- 
ence vs. undersampled data. Linear regression was per- 
formed, resulting in a slope of 1.01 (±0.04 95% 
confidence interval) and an offset of -2.01% of the range 
of contraction (±4.79% 95% confidence interval). 

In-vlvo measurements 

Figure 5 shows short axis slices reconstructed from fully 
sampled and 3- and 4-fold retrospectively and prospect- 
ively undersampled data. Frames at 27 ms after detection 




■R3BR4nR5nR8 

Figure 4 Comparison of circumferential shortening and 
rotation for fully sampled reference data and retrospectively 
undersampled data from computer simulation. Relative 
differences are presented as average over tine entire LV (a-b) along 
with slope and offset of linear regression and the corresponding 
correlation coefficient and the standard error of the estimate 
(SEE) (c-j). 

V J 




Figure 5 Comparison of short-axis views reconstructed from 

fully sampled (ref) and 3- and 4-fold retrospectively 

(retrospect) and prospectively (prospect) undersampled data. 

Three different time points are shown. 
V . J 



of the R-wave, at end systole (279 ms) and at mid dia- 
stole (559 ms) are shown. 

Linear regression analysis performed on fully sampled 
and retrospectively undersampled data is shown in 
Figure 6. For each regression the equations of the linear 
fit, the correlation coefficients and the SEE are given. 
Figure 6a-d show the analysis of polar coordinates of 
tracked points and Figure 6e-h demonstrate sector-wise 
comparison of time curves of circumferential shortening 
and rotation for R = 3 and 4. 

In Figure 7 shows the Bland- Altman comparison of 
fully sampled and retrospectively undersampled in-vivo 
data for R = 3 and 4. Dashed lines represent the mean 
difference (light grey) and the 95% limit of agreement 
(2SD) (black). In Figure 7a-d provides the comparison of 
radial coordinates (r = radius and (|) = angle) of tracked 
contour points. Mean differences ± 2SD of r were 
0.0 ± 3.3% and -0.2 ± 4.1% for R = 3 and 4 and -0.2 ± 1.9° 
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Figure 6 Linear regression analysis for point-wise and sector-wise comparison of fully sampled (ref) and retrospectively undersampled 
(und) data in healthy volunteers (R = 3 and 4). Point-wise comparison was performed for polar coordinates r (a, c) and cp (b, d). For 

circumferential shortening (e, g) and rotation (f, h) entire motion curves were compared per sector across the entire LV. 



and -0.3 ± 2,2° for (|). In Figure 7e-h Bland- Altman plots 
for peak circumferential shortening and peak rotation 
are given. The mean differences ± 2SD for circumferen- 
tial shortening were -0.2 ± 4.1% and -0.1 ± 4.2% for R = 3 
and 4. For rotation, mean differences ± 2SD were 
0.5 ± 1.8° and 0.7 ± 1.7° for R = 3 and 4. The mean differ- 
ences of torsion were 0.45 ± 2.22° (p = N.S.) and 
0.05 ± 2.24° (p = N.S.) for R = 3 and 4, respectively. 

Figure 8 shows the resulting Bland- Altman analysis of 
peak motion (a-d) and correlation between full time 
curves of motion (e-h) comparing fully sampled data 
with data acquired with prospective undersampling 
(R = 3 and 4). Mean differences ± 2SD of peak circumfer- 
ential shortening were -0.1 ±5.2% and -0.7 ±6.2% for 
R = 3 and 4. Mean differences of peak rotation were 
found to be 0.44 ± 1.80° and 0.73 ± 1.67° for R = 3 and 4 



while differences in torsion were 0.48 ± 4.20° (p = N.S.) 
and 0.03 ± 4.48° (p = N.S.) for R = 3 and 4, respectively. 

Bulls -eye plots of peak circumferential shortening and 
profiles of circumferential shortening and rotation 
derived from undersampled (R = 4) and fully sampled 
data are compared in Figure 9 for one healthy subject. 
Spatially depend differences are not observed. 

Figure 10 compares circumferential shortening 
obtained from undersampled and fully sampled refer- 
ence data for the five patients in relation the LGE 
findings. Differences in peak circumferential shortening 
in sectors with more than 50% of their area present- 
ing LGE and sectors with less than 50% LGE were 
statistically significant in all patients. Mean differ- 
ence ±2SD in torsion between undersampled and fully 
sampled reference data in patients was 0.72 ±2.14°. 
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Figure 7 Bland Altman plots of point-wise and sector-wise comparison of fully sampled and retrospectively undersampled in-vivo data 

(R = 3 and 4) in healthy volunteers. Mean differences (grey) and 95% levels of agreement (2SD) (black) are indicated by the dashed lines. 
Analysis was performed on the entire LV. Point-wise comparison was performed for the polar coordinates r (a, c) and cp (b, d) for each tracked 
point. Circumferential shortening (e, g) and rotation {f, h) were compared per sector. 
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Figure 8 In-vivo results of the comparison between fully sampled reference data (ref) and data acquired with an acceleration rate of 
R = 3 and 4 (und). Mean differences (grey) and 95% level of agreement (2SD) (black) are shown for peak circumferential shortening (a, c) and 
rotation (b, d) of the entire left ventricle. Motion curves of circumferential shortening (e, g) and rotation (f, h) of undersampled and reference 
data are plotted against each other and linear regression was performed. 



The SNR in the fully sampled in-vivo data was 
36 ± 12 on average. 

Discussion 

In this work undersampled 3D CSPAMM tagging in 
combination with k-t PCA has been implemented and 
validated. 



The analysis of simulation results has revealed a max- 
imum applicable acceleration factor of four. Generally, 
good correlation was found between fully sampled and 
retrospectively undersampled data over the entire left 
ventricle. Analysis of relative differences in circumferen- 
tial shortening showed a slight decrease in error going 
from the apex towards the base of the heart. As motion 




Figure 9 Comparison of circumferential shortening for fully sampled data and 4-fold undersampled data. The first column shows maps of 
peak circumferential shortening for fully sampled b) and undersampled data c). Resulting motion profiles taken from an equatorial slice are shown 
for circumferential shortening (d, e) and rotation (f, g). The dots represent the actual data points while the line represents the polynomial fit. 
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Figure 10 Maps of peak circumferential shortening (csh) for five patients with myocardial infarction. The top row represents results from 
fully sampled data while the second row shows data from prospectively undersampled acquisitions. The third row shows the result of late 
gadolinium enhancement imaging (LGE). The myocardial signal intensity post contrast was averaged over each sector and normalized to the 
signal intensity from completely infarcted tissue. The bottom row shows box-plots of peak circumferential shortening for sectors with more than 
50% late gadolinium enhanced signal (grey) and sectors with less than 50% late gadolinium enhanced signal (black). The box represents 50% of 
the data points and the error bars 90%. Statistical significance is indicated by # for p < 0.002 and * for p < 0.0002. 



was normalized for each slice, this observation is asso- 
ciated with the low magnitude of contraction of apical 
myocardium. 

Different transmural extensions of infarcted tissue 
were investigated by computer simulation. Motion ana- 
lysis showed a good correlation between undersampled 
and fully sampled reference data. Transmural reduction 
of contraction as low as 20% could be distinguished 
from healthy fully contracting tissue. Although fully 
infarcted tissue was simulated as well, sector-wise ana- 
lysis always showed contraction greater than 6%. This is 
due to the choice of the position and size of the sectors 
used for analysis. Sectors with non-contracting tissue 
contained slightly contracting neighbouring tissue, which 
leads to a non-zero average circumferential shortening 
within a sector. 

In vivo comparison between fully sampled data and 
retrospectively undersampled data showed good agree- 
ment as the mean difference was less than 1% (% of 
range of motion). Differences between reference data 
and prospectively undersampled acquisition were found 
to be larger in comparison to differences seen relative to 
retrospectively undersampled data. This finding is 



related to multiple issues. On one hand, reference data 
and undersampled data were acquired in two consecu- 
tive scans. Despite the use of a gating respiratory naviga- 
tor, different breath hold levels are possible and slices 
used for analysis might hence differ in position. On the 
other hand the fully sampled reference data itself was 
acquired in three consecutive breath holds, which, des- 
pite the use of a respiratory gating window of 5 mm, can 
lead to stack misalignment within the reference data. 
The resulting circumferential shortening is dependent 
on the transmural position of the tracked contour. If 
stack misalignment occurs a contour that appears to be 
midmural in one stack can be positioned more epicar- 
dially or endocardially in the consecutive two stacks 
from which the motion orthogonal to the first stack is 
derived from. Hence a contour point is tracked at three 
different transmural positions. Therefore the motion 
profiles obtained from fully sampled data, used as gold 
standard reference in this work, might have been com- 
promised. Image registration was not performed in this 
study. The implementation of tailored registration algo- 
rithms that can handle orthogonal tagging modulations 
requires further investigation. 
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Maps of peak circumferential shortening as well as the 
corresponding motion profiles were found to agree well 
for reference and undersampled data. The magnitude of 
circumferential shortening is comparable to previously 
reported values [24,25,50,51]. Besides potential motion 
smoothing, expressed by the positive mean difference 
for rotation, a slice by slice comparison of prospectively 
undersampled and fully sampled data might have been 
biased by an offset in breath hold level for both 
acquisitions. 

In this study, standard 2D single slice tagging data 
were not available for reference. In order to capture 
through-plane motion and hence make data comparable 
to 3D tagging, non-standard extensions of 2D tagging 
such as slice following [15] or zHARP [18] are required. 
Accordingly, assessment of error in the present work 
was relative to fully sampled 3D CSPAMM data only, 
which is a clear limitation. Nevertheless, there have been 
previous records of validating 3D CSPAMM on healthy 
subjects and patients, which may serve as benchmark 
[24,25]. 

Spatial resolution in each of the three orthogonal 
stacks in 3D CSPAMM may be considered coarse 
(3.8 X 5.4-7.2 X 7.2 mm^). The resolution given in readout 
direction is, however, directly linked to the tagline spa- 
cing when using HARP analysis. Higher temporal reso- 
lution (< 30 ms) is desirable as this allows separating 
data in x-pc space further and hence improves recon- 
struction accuracy in k-t PC A. However, this results in 
more heart phases and hence more RF excitations per 
cardiac cycle reducing the contrast-to-noise ratio of the 
tagged data. 

Five patients with myocardial infarction were exam- 
ined using undersampled tagging and LGE imaging. 
Maps of circumferential shortening derived from refer- 
ence data and prospectively undersampled data agreed 
well and infarcted regions could be localized. Direct cor- 
relation of peak circumferential shortening with the area 
of delayed enhancement was found to be only moderate. 
A potential limitation of infarct detection is related to 
averaging of motion within sectors. A six-sector per slice 
model [24,25] results in 16.7% of total myocardial mass 
captured per sector and per slice. For example, non- 
viable tissue in border zones of infarction is passively 
moved and compressed and hence the extent of motion 
abnormality may be overestimated. This issue may be 
addressed by increasing the number of sectors per slice 
and by using multiple circumferential contours covering 
the entire transmural extent in future work. 

Several strategies may be envisioned to shorten the 
relatively long breath hold durations (-20 sec) of the 
present protocols. First, sampling of training data can be 
fully integrated into the acquisition of undersampled 
data. Such a variable density EPI approach provides 



shorter overall scan duration at the expense of reduced 
temporal resolution. Second, the separation of training 
data and undersampled data allows splitting the data ac- 
quisition into two breath holds. Compared to the 
method proposed by Rutz et al. [24] the three stacks are 
not acquired within three consecutive breath holds, but 
could be acquired within one breath hold for high reso- 
lution undersampled and a second breath hold for low 
resolution training data. Since the training data has very 
low spatial resolution in the phase-encode directions 
(2.3 X 2.3 cm^), differences in breath hold levels between 
acquisitions no longer pose a concern. 

Conclusion 

Undersampled cine 3D tagging in conjunction with k-t 
PCA reconstruction enables significant reduction in scan 
time of whole-heart tagging and facilitates efficient and ac- 
curate quantification of shortening, rotation and torsion 
of the left ventricle. Using 3-fold undersampling the entire 
cine 3D tagging acquisition could be accommodated in a 
single breath hold and feasibility in volunteers and 
patients was demonstrated. Future work is dedicated to 
shorten breath hold durations further and to apply the 
method in larger patient cohorts to prove clinical value. 
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